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A quantization procedure without Hamiltonian is reported which starts from a statistical ensemble of 
particles of mass m and an associated continuity equation. The basic variables of this theory are a 
probability density p, and a scalar field S which defines a probability current j — pVS/m. A first equation 
|' for p and S is given by the continuity equation. We further assume that this system may be described 

by a linear differential equation for a complex-valued state variable X- Using these assumptions and the 
simplest possible Ansatz x(p, S), for the relation between \ and p, S, Schrodinger's equation for a particle 
of mass m in a mechanical potential V(q, t) is deduced. For simplicity the calculations are performed for 
' a single spatial dimension (variable q). Using a second Ansatz x(p, S, q, t), which allows for an explicit 

q, i-dependence of x> one obtains a generalized Schrodinger equation with an unusual external influence 
described by a time-dependent Planck constant. All other modifications of Schrodinger' equation obtained 
within this Ansatz may be eliminated by means of a gauge transformation. Thus, this second Ansatz may 
be considered as a generalized gauging procedure. Finally, making a third Ansatz, which allows for a non- 
unique external q, i-dependence of x, one obtains Schrodinger's equation with electrodynamic potentials 
A, (f> in the familiar gauge coupling form. This derivation shows a deep connection between non- uniqueness, 
quantum mechanics and the form of the gauge coupling. A possible source of the non-uniqueness is pointed 
out. 



P.A.C.S.: 03.65.Ta; 06.20.Jr; 11.15.-q 
1 Introduction 

Usually a physical system, which one wants to describe quantum mechanically, is first identified in the 
context of classical physics and then somehow transferred to the quantum mechanical domain; this 
process is referred to as " quantization" . Often the first step in the quantization process is the tacit 
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assumption that a Hilbert space is associated with the examined system. Then, the remaining task is to 
find the proper algebra of operators. A more direct method which avoids this assumption, is to "derive" 
Schrodinger's equation, i.e. to find premises which imply Schrodinger's equation. This may be done in 
several ways. The method which was historically at the beginning of quantum mechanics [21] ("wave 
mechanics" ) starts from the Hamilton- Jacobi equation of a classical system and tries to deduce from it 
- with the help of suitable modifications [31 21 [7] - the corresponding quantum-mechanical equation for 
the time development of the system. Other premises leading to Schrodinger's equation include special 
assumptions about the structure of momentum fluctuations [5J and the principle of minimum Fisher 
information [15] . 

In this paper, a new quantization procedure is reported which shares with the last two examples the 
property that it does not start from a single-particle picture but from a statistical ensemble. The simplest 
nontrivial system, a spinless particle of mass m in nonrelativistic approximation is investigated. In order 
to define the subject of this work more precisely we start from the classical Hamilton- Jacobi equation for 
the action function S(q,t), which depends on the particle coordinates qk and the time t. It is given by 



dS(q,t) | 1 S dS(q,t) 
dt 2m \ dq 



V(q,t) = 0, (1) 
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if the movement takes place under the influence of a potential V{q, t). The momentum field, that appears 
in Eq. ([T]) is given by 

Pk m = 9 -f^. (2) 

oqk 

The fact that the Hamilton-Jacobi equation is the ideal starting point for the transition from classical 
physics to quantum mechanics is formally based on the following well-known reformulation of the time- 
dependent Schrodinger equation 

~m t) + V(q, t)m t) = ^^(q, t). (3) 
If the complex-valued variable ip(q,t) is, without any restrictions of generality, written in the form 

mt) = ^/pW)^ sm \ (4) 

then one obtains from Eq. (JU) , by calculating the real parts of both sides, the relation 

dp(q, t) dp dS(q, t) 



dt dq m 



0, (5) 



which is a classical (no h occurs) continuity equation for the probability density p and the probability 
current pp/m. Equating the imaginary parts of both sides of Eq. ([3]) one obtains the relation 



^4»)- +f(i „4« (6) 

dt 2m \ dq J 2m y/p(q,t) 

which differs from the Hamilton-Jacobi equation ([I]) only by the single term on the right hand side. 
Eq. ([(j]) is sometimes referred to as quantum Hamilton-Jacobi equation. 

Due to this similarity there have been attempts [H [7] to use Eq. ([1]) as a starting point and to derive 
Eq. ([3]), which presents the basis of quantum mechanics, by justifying introduction of the crucial quantum 
term appearing in Eq. ([6|). In the present work we go a different route, starting from assumptions which 
are simpler in certain respects. We postulate the existence of a statistical ensemble but do not start from 
the Hamilton-Jacobi equation itself. Instead, our basic postulate is the validity of a continuity equation 
(of the above type) , interpreted as a local conservation law of probability. Since we have two unknown 
functions (p und S) and only one single equation, we clearly need further assumptions - and a second 
differential equation - in order to arrive at a mathematically well-defined problem. Our second assumption 
is very simple and of a purely formal nature. We require that both equations - the one already known 
and the second one still to be found - may be expressed mathematically as a single equation for a single 
complex state variable x- This second assumption expresses something like the postulate of maximal 
mathematical simplicity. As we know, this postulate may be quite successful in physics, in particular if 
combined with other ideas. 

Thus, the continuity equation has to be "extended" to the complex domain. This task may be 
described briefly as follows: Consider a complex-valued variable x which depends in an unspecified way 
on the real variables S and p. Which differential equations for x exist, whose real (or imaginary) part 
agrees with the continuity equation and which functional dependencies x(/9, S) are compatible with this 
requirement ? Note that both the functional dependence of x an d the shape of the differential equation 
are unknown "variables" of this problem. A detailed formulation of this problem is reported in the next 
section [2]). The calculation, reported in section [3]) and appendix lAl has been performed for simplicity 
for a single spatial dimension and for a reduced class of differential equations obeying several additional 
constraints. The result is Schrodinger's equation for a particle in an external mechanical potential. 
Further, we formulate and justify in section [3|) the conjecture that all additional constraints except 
linearity may be omitted. 
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In section 0]) the original Ansatz is extended by allowing for an additional, explicit space-time de- 
pendence of the state variable x- This is our first attempt to derive the minimal coupling rule, which 
is obeyed by essentially all fundamental interactions, in the present context. It turns out (details of 
the calculation are reported in appendix [B| that a second, very unusual external influence, besides the 
potential V, appears in the Schrodinger equation. It takes the form of a time-dependent Planck constant. 
All other modifications due to the extended Ansatz are spurious, because they may be eliminated from 
the Schrodinger equation by a gauge transformation. The gauge field itself cannot be derived by means 
of this Ansatz. In section [5]) our second attempt is undertaken to derive a gauge field. The Ansatz of the 
last section is once more extended by allowing for a non-unique space-time dependence of x- The result 
is Schrodinger's equation for a charged particle in an external electromagnetic field. 

Section [6]) contains a detailed discussion of all assumptions and results and may be consulted in a first 
reading to obtain an overview of this work. It also contains remarks of a speculative nature concerning 
the relation between the classical theory of charged particles and fields on the one hand and the form of 
quantum mechanics and gauge coupling on the other hand. In the last section [7]) one finds concluding 
remarks. 

2 Formulation 

We consider a statistical ensemble which is described by a probability density p(q, t) and a probability 
current j — pp/m. We assume that the momentum field p may be written as the gradient of a scalar 
function S(q,t), as in this simplest possible form of p holds in Hamilton- Jacobi theory [20 . Our 
ensemble can thus be described by two variables, the real fields p and S. From the continuity equation 
we obtain immediately Eq. (|5|) as a first differential equation for p and S. The physical meaning of the 
function S(q,t) is still unclear. Of course, in the classical limit of the theory to be constructed, S(q,t) 
should agree whith the action function of classical mechanics. 

We are confronted with the problem of finding a second differential equation for p and S. Two 
solutions of this problem may be found in nature, classical statistical mechanics and quantum theory. 
Let us compare these two solutions trying to learn something about the quantization process. The second 
equation is {1} in classical physics and ([6]) in quantum theory. One sees immediately that in quantum 
theory the additional term in (J6j> implies a dramatic change of the basic concepts. This term is not 
an externally controlled input parameter but describes a coupling between the identity of the examined 
object and the statistics. But let us put aside these physical aspects. Instead let us ask the more formal 
question in which of the two theories the task of formulating the basic equations has been solved in a 
simpler way. This is certainly quantum theory, because the two differential equations ([5|) and (J6j> can 
be combined to a single equation, the Schrodinger equation (|3|) which has on top of that a simple linear 
structure. A similar unification is not possible in classical physics (one finds [20) a nonlinear equation 
containing both ip and ip*)- This situation suggests that this principle of simplicity may be used, turning 
the logical direction around, to derive quantum mechanics from classical mechanics or to use it at least 
as an essential part of the quantization process. 

This principle of simplicity implies that the two differential equations for p and S, namely Eq. (|5|) and 
the second one which is still to be found, may be combined into a single equation for a two-component 
variable Xi which we assume to be a complex quantity. The variable x replaces p and S and should of 
course be a function of p and S. If we define x hi terms of its real and imaginary parts (which has the 
advantage of linearity in comparison with the polar representation) x takes in the simplest case (more 
general relations will be considered in the next sections) the form 

X(P,S) = xi(p,S) +iX2{p,S). (7) 

We assume furthermore that the differential equation for Xi which will be referred to as "state equation" , 
is a partial differential equation (with complex coefficients). The independent variables of this equation 
must be the usual space-time variables qk, t, which also occur in the continuity equation. 

Each differential equation for x whose real (or imaginary) part agrees with the continuity equation 
represents a possible extension of the latter. Its imaginary (or real) part provides us automatically with 
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a second differential equation for our two variables. Thus, each one of these equations defines, in a purely 
formal sense, a physical theory motivated by the principle of simplicity. We denote the class of differential 
equations defined in this way by C. We know that C is not empty, ft contains certainly Schrodinger's 
equation ([3]); it is a calculation of few lines to derive the continuity equation (O starting from ([3]). But 
the inverse problem, to start from the continuity equation {5J and to derive the set of all extensions, i.e. 
the set C, may be less trivial. In particular one cannot assume that C contains only a single equation. 

Being interested in the transition to quantum mechanics, we want to know under which conditions 
Schrodinger's equation can be derived. For each quantization procedure the important question to ask is 
how the final result is obtained. For that reason we make no use of symmetry considerations of any kind 
in this paper. Thus, which additional conditions are necessary in order to single out, from the set C, 
Schrodinger's equation as only remaining equation ? We introduce for brevity a symbol, say A, to denote 
the set of all additional conditions defined in this way. C and A define together a quantization procedure. 
Obviously, this quantization procedure will only be convincing for a small number of (physically appealing) 
additional conditions, i. e. if A is "small" . What we are looking for is the smallest set of such assumptions. 

If A should turn out to be empty, this would result in a very impressive quantization procedure; 
it would mean that quantum mechanics could be derived from only two assumptions, the continuity 
equation and the existence of a complex state variable [as well as the tacit assumption that the laws of 
nature may be formulated as differential equations of the conventional type, compare section [6j] . In the 
next section it will be shown that this is not the case. However, the results indicate that only a single 
additional condition, the linearity of the differential equation, is required. 

How to find concretely the set A ? This set A is the smallest set of assumptions which eliminates 
all equations except the Schrodinger equation from the set C. An obvious strategy is to start from 
a set of strong assumptions, say Ai, which defines a small class of differential equations (all beeing 
loosely speaking "similar" to Schrodinger's equation). If the conditions implied by Ai and C single 
out Schrodinger's equation, one may eliminate some of the conditions contained in Ai and test if the 
corresponding, smaller set A2 of conditions may be used instead to lead to the same result. The final 
solution is the smallest possible set A obtained in this way. 

In this paper a complete solution to this problem is not given. All explicit calculations are performed 
using a set Ai, defined by the following constraints for the state equation: 

• linearity, 

• a single space dimension (variable q), 

• only derivatives of first order in t , 

• only derivatives up to the second order in q, 

• no mixed derivatives. 

However, the structure of the solutions leads to the conjecture on A mentioned above, which is discussed 
in more detail in the next section. 

3 No interaction 

In the framework of the strategy set up in the last section the possibility of interaction, i. e. of an 
external influence on our single particle system, was not taken into account. In principle there are two 
possibilities to do that; we may either modify the set A or the definition of \- In this work the second, 
more general, possibility will be chosen. In this sense, the simplest Ansatz x = x{Pi S), dealt with in the 
present section, describes a situation without interaction. The meaning of the term interaction defined in 
this way may differ (and does in fact differ) from the conventional one; we continue to use it for simplicity. 

We now have the following concrete mathematical problem: Find all differential equations obeying 
the conditions Ai whose real part (or imaginary part - this can be fixed arbitrarily and does not represent 
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a real constraint) agrees - after splitting off an arbitrary multiplicative factor - with the one-dimensional 
continuity equation. Our unknown variables are the coefficients of the differential equation, the state 
function x, an d the factor F. 

The basic complex- valued dynamic variable x(p> S) is a function of the real variables p(q, t) and S(q, t) 
and is written in the above form ([7]). The multiplicative factor may be an arbitrary complex function 
F(p, S, q, t). The coefficients of the linear differential equation are denoted by a, b, d, e and are arbitrary 
complex functions of q, t. Then, the fundamental requirement implied by C and Ai takes the form 
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1 ^ a wt x+h h x+d h x+ex 



dt dqm dq m dq 2 



The right hand side of ([5]) is given by the continuity equation. 

We have only a single equation |8j and many unknowns, namely the complex functions S), a(q, t), b(q, t), d(q, t), e(q, t) 
However, Eq. ([5]) is a very strong requirement because the variables p, S are arbitrary. Consequently, the 
coefficients of p, S - and of all derivatives of p, S - must agree on both sides of ([8]). Using this fact and 
introducing real and imaginary parts of a, . . . , e and F according to a = a± + ta 2 , ■ ■ ■ , e = e\ + ie 2 and 
F = Fi + iF 2l we obtain from Eq. ([5]) the following 10 conditions 

(a 1 F 1 - a 2 F 2 )^ - (a 2 F 1 + ai F 2 )^=0, (9) 

(diFi - d 2 F 2 )-^- - {d 2 F x + dxF 2 )-^=0, (10) 

(diFi - d 2 F 2 )^- - (d 2 Ft + d!F 2 )-^-=0, (11) 
(biF x - b 2 F 2 )^- - (6 2 iq + b 1 F 2 )^-=Q, (12) 

(fuFi - hF 2 )^ - {b 2 Fx + hF 2 )^=0, (13) 
(diFi - d 2 F 2 )^ - {d 2 F x + d 1 F 2 )^-=Q, (14) 

(oiFi - a 2 F 2 )^- - (a 2 Ft + a x F 2 )^=l, (15) 
op op 

(*'--*»)^-<*« + *«»^-s:' (16) 

(<ZlF, - <i 2 ft)% - (*fi + (17) 
ob Ob m 

(eiFi - e 2 F 2 ) X i - (e 2 iq + ei F 2 ) X2 =0. (18) 

We have 12 unknown quantities and 10 equations, but the unknown variables a±,a 2 , . . . ,e±,e 2 and F\, F 2 
occur in ©- (dU) only in the combinations 

ai = aiFi — a 2 F 2 , a 2 = a 2 Fi + aiF 2 , (19) 

h=b l F l -b 2 F 2 , b 2 = b 2 F 1 +b 1 F 2 , (20) 

dj = d 1 F 1 - d 2 F 2 , d 2 = d 2 Fi + d x F 2 , (21) 

ei = ei F! - e 2 F 2 , e 2 = e 2 F x + e x F 2 . (22) 

Thus, we have 10 equations for 10 unknown quantities ai, a 2 , . . . , ei, e 2 , XiiX2- Most of the equations (J9j>— 
(fig)) look like differential equations. However, the coefficients apearing in these equations do also belong 
to our set of unknown variables. Thus, relations ([9])- (fT8|) may either be used as differential equations for 
the unknown functions \i (P; S) > X2 (p, S) ... or - if these functions are already known - as constraints for 
the unknown coefficients. In the latter case, one has relations which must hold for arbitrary p, S, i.e. for 
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each one of the considered conditions [out of the set ©- (JTSJ)] the coefficients of all linear independent 
functions of p, S must vanish separately. While this type of problem may be somewhat unusual, its 
solution is straightforward and does not require any unusual mathematical methods (its solution could 
possibly be simplified by using more sophisticated methods). 

The first part of the calculation, the determination of \ii Fi is described in detail in appendix lAl The 
result is 
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d\ cos 
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C2 

2m\df 
;S + C 5 



;S + C 5 



rS + C 6 



rS + a 



-S+C 6 



(25) 



(26) 



where c±, c 2 are linear combinations of coefficients, defined in Eq. (|114p . and \d\ is the modulus of the 
complex number d. 

Only a part of all conditions, namely ([9]), (fTOj) . (fi"4"]) . ([15]), (fl6)) . and (fl7|) have been used in the course 
of the calculations leading to Eqs ([23ll - (|26p (see appendix IX| . In these conditions only the coefficients 
a, and di appear. Conditions (|12p . ()13|) . containing the 6i, and condition (|18|) . containing the e%, have 
not been used. The conditions not yet used play the role of constraints for the constants of integration 
C3, C4, C5, Cq and the real- and imaginary parts of the coefficients a, b, d, e. We will, for brevity, refer 
to the totality of all these quantities as "parameters" . 

Our 12 parameters are constants with respect to the variables p, S, but may be arbitrary functions of 
q, t, as far as the calculation reported in appendix [Al is concerned. However, they have to obey the basic 
requirement of this section, that \ depends only on p, S, and not on q, t. This implies that the following 
parameters or combinations of parameters 



C3, C4, C5, Cq, 



('2 



2m\d\ 2 ' 2m\d\ 2 



(27) 



are constants (we use this term now in the usual sense of being independent of q, t). 

In order to obtain the explicit form of the constraints for the parameters, Eqs. (|23 |) - ([26l) have to be 
inserted in the conditions not yet used ([TTJ) , (fl2^) , (fT3"j) and (TTg|) . These conditions hold for arbitrary values 
of p, S. This leads to C3 = C4 = while no constraints exist for C5, Cq. For the coefficients one obtains 
the conditions 



ci = aidi + a 2 d 2 = 0, 

dibi + d 2 b 2 = 0, 

d\b 2 - d 2 b\ = 0, 

a 2 di - a\d 2 = r% (d\ + d 2 ) 

e 2 di - eid 2 = 0, 



(28) 
(29) 
(30) 
(31) 
(32) 
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where r\ in (f5T|) is an arbitrary real constant. As a consequence of (|2"5)) - (n?2l all parameters (with the 
exception of the arbitrary real constants C5 und C§ ) may be expressed in terms of three arbitrary real 
functions di(q,t), ^2(5, t), f(q,t) and the real number r±. In terms of the complex coefficients one 
obtains the following result: d is an arbitrary complex-valued function of q und t, b = 0, and a and e are 
determined by d according to the linear relation 

a = vr\d, e = fd. (33) 

Now, all conditions have been taken into account. The complex quantities \ an d F may we written as 

x= _^ e - C8 Vpe'(& s+Cs ), (34) 
n 

F=l }H e c^{^s + c,) (35) 
md 

and the differential equation for \ takes the form 



•:'4)- • </0 • (•%) 



V(q,t)1>, (38) 



If S 1 obeys the continuity equation ([5]), then it has the dimension of an action. Consequently, the 
constant 2m/r\ has the dimension of an action; we identify this constant with Planck's constant h. The 
constant d may be canceled and the constant e has the dimension cmT 2 . We may now come back to a 
more conventional notation by introducing quantities tp(q, t) and V(q, t), which replace x(<li t) and f(q, t) 
and are defined by 

^^V^Pe*^, V(q,t) = - — f(q,t). (37) 
Then, Eq. ([36]) takes the form 

7i(?^> h 2 d 2 ifj 
1 dt 2m dq 2 

which agrees with the one-dimensional Schrodinger equation for a particle of mass m in an external, 
mechanical potential V(q,t) (the terms containing C5, Cq, which correspond to the usual freedom in 
amplitude and phase, have not been written down). 

Thus, using the Ansatz \ — x{pi S) an< ^ the set of conditions Ai we found that only the one- 
dimensional Schrodinger equation yields the continuity equation as its real part. More precisely, we 
found an infinite number of such equations differing from each other by the choice of an arbitrary func- 
tion V(q, t). We obtained a coupling to an external potential but no coupling to a gauge field; the 
meaning of our term "no interaction" should be adjusted accordingly. If the interaction is considered as 
a secondary aspect, the above assumptions define a quantization procedure, i.e. a way to perform the 
transition from classical physics to quantum mechanics. 

Can we improve this quantization procedure by reducing the number of assumptions contained in Ai ? 
Let us keep first the linear structure as well as the one-dimensionality of the equation and allow for higher 
derivatives with respect to t und q. Thus, new terms appear in the fundamental condition l|8|). Will these 
new terms survive or will the corresponding coefficients vanish? All derivatives of p and S, which occur 
in the continuity equation, on the right hand side of ([5]), appeared already in the differential equation 
defined by Ai. Therefore, all the new terms, due to the higher derivatives, will have to vanish (the 
corresponding coefficients are all solutions of homoneneous equations) just as, in the above calculation, 
the term proportional to the first derivative with respect to q had to vanish. 

As a next step, let us allow for three spatial dimensions and keep only the condition of linearity. To 
get an idea what happens in this case, preliminary calculations using, instead of ((8]), the fundamental 
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condition 




) 



(39) 



have been performed (summation from 1 to 3 over double indices). The differential equation displayed 
in |39|) contains all derivatives of qu,t up to second order taking into acount also the possibility of 
anisotropic coefficients. These (not completed) calculations indicate (i) that the anisotropy of the coeffi- 
cients is, as one would expect intuitively, suppressed by the isotropy of the continuity equation, (ii) that 
the coefficients of the higher derivatives vanish, for the reason just mentioned, and (iii) that an exact calu- 
lation starting from (|39[) . which is completely analogous to the one above, leads to the three-dimensional 
Schrodinger equation (J3]). 

The remaining condition in Ai is the linearity of the differential equation. It can not be eliminated; 
one realizes immediately that e.g. a cubic nonlinearity of the form r\%p\ 2 ip, where r is real, is compatible 
with the continuity equation. Thus we conjecture that our final quantization procedure is defined by a 
smallest set A, which contains only a single condition, namely the linearity of the differential equation. 

What is the physical meaning of the remaining condition of linearity ? A plausible explanation is that 
linearity is a consequence of the probabilistic interpretation [9J of the wave function, which excludes (in 
general) predictions about single events. There is a large literature on the relation between linearity and 
indeterminism; we only mention Mackey's axioms for statistical theories [14] and Caticha's derivation 
of quantum mechanics from the rules of manipulating probability amplitudes pQ. On the other hand 
it is well-known that macroscopic quantum phenomena exist in nature, such as superconductivity and 
superfluidity, which may be successfully described by nonlinear terms in the corresponding complex state 
variables. This is not in conflict with the above explanation. In contrast to the true wave function, the 
state variables of these many-body theories do not possess a probabilistic interpretation. They lost their 
"immaterial" (probabilistic) meaning in the thermodynamic limit and may be directly measured in single 
experiments. The condition of linearity is required to exclude such theories from the consideration. 

4 First attempt to derive interaction 

The Ansatz \ — x{Pi S) used in the last section led, somewhat unexpectedly, to a term Vip in 
Schrodinger's equation (|38p . describing an interaction by means of an external mechanical Potential 
V. A term describing coupling to a gauge field did, not unexpectedly, not appear. In this section we 
start to study the following question: Is it possible to obtain this type of coupling, which is obeyed by 
all fundamental interactions, in the present framework ? If possible it requires, at any rate, a different, 
properly generalized Ansatz. 

There is an obvious possibility to generalize our Ansatz with regard to an external influence: we may 
allow for an additional explicit g, ^-dependence of our state variable \i i- e - write 



instead of \ — x(Pi S). Such an extension seems quite natural if one wants to describe an external, 
otherwise unspecified influence on the system described by \- We will study the consequences of (|4T))) 
keeping all other assumptions unchanged. 

Using (|40|) does not change the appearance of the basic condition ([8|) in a fundamental way but 
additional derivatives with respect to the explicit q, ^-dependence have to be taken into account. To 
indicate this difference, we replace the symbols for the partial derivatives by symbols for total derivatives. 
Thus, our fundamental condition, generalized according to the new Ansatz (|40[) takes the form 



X = X(P, S, q, t), 



(40) 




(41) 
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The form of the factor F(p, S, q, t) and the coefficients a, 6, d, e remain unchanged; the latter may be 
arbitrary functions of q,t. The problem defined by (|4i[) contains now, in comparison to ([5]), an even 
larger number of unknown functions. Fortunately, it will turn out that the equations to determine the 
p, S'-dependence and those for the q, i-dependence are "decoupled" in the sense that they may be solved 
one after the other. 

Comparing the coefficients of the derivatives on both sides of Eq. (|4"Tj) we now obtain the following 
fO conditions 



_ dxi _ dx2 „ 



(42) 



di?£-dfe=0, (43) 



dp 2 2 dp 2 
d 2 Xl - d 2 X2 _ n 

d X i - 8x2 . /t d 2 Xi t d 2 X2\ „ 

& 2^— + 2 [di—- — d 2T — = 0, 

Op op V Opoq Opoql 

,- %i ,- %2 , „/ T d 2 xi 7 d 2 X2 \ n 
~ + 2 { dl dSo- q - d2 dSd- q ) = °- 



(44) 
(45) 
(46) 

- dxi 7 dX2 n , , - , 

ai— "2^— =0, (47) 



"2^— 

ap op 

_ 9xi _ 9x2 . 

a 2^ — =1 > 

op ap 

r 9 2 xi _ t 9 2 X2 _ 1 

'dpds 2 d P ds 2m' 

7 9xi 7 9X2 _ p 

. -9xi _ 9x2 . ,- 9xi ,- 9x2 . 7 9 2 xi r 9 2 X2 „ 

eiXi - e 2X2 + a i^rr - a 2~z~r + h — b 2 -^ h d 1 — T - - c^^r^- = 0. 

Ot ot oq Oq oq 2 Oq 2 



(48) 
(49) 
(50) 
(51) 

Again, each one of the conditions (|42^ - (f5Tj) may split into several sub-conditions, if on the left hand 
side several linear independent functions of p, S occur. Comparing with the previous "interaction-free" 
conditions (|9"j)- (|18[) . we see that seven of the Eqs. (|42 p - (|?T|) agree with corresponding equations in the 
set ©-([IH]). Only gSJ, (gHD, and differ from the corresponding previous conditions (|f 2[) . (jf 3|) . 
and (|f 8p by new terms. These new terms appear as a consequence of the explicit q, ^-dependence of Xi 
and lead obviously to a coupling of previously independent coefficients. 

A great simplification of the problem defined by ([42 |) -([5T |l takes place as a consequence of the fact that 
conditions (fT2|) . (fl~3| and (fl8|) have not been used during the calculation of the p, S'-dependence of x hi 
section (only conditions containing a and d were required). Therefore all conditions in the set (jl2")) - (f5Tj) 
which are necessary to obtain this relationship remain unchanged. Since the additional q, i-dependence 
does not affect this part of the calculation, the formal results for x and F from section [3]) may be taken 
over without modification for the present calculation. Thus, for the real and imaginary parts of x we 
obtain 

Xi (p, S, q, i)=-^!^cos ( -^S + C 5 ) e -^ S - Ce + C 3 , (52) 



c 2 V 2m M 5 

2m\d\ 2 ( c 2 



X2 (p,S,q,t)= ^VP^[2^S + C 5 )e ^ + C 4 . (53) 



r S — Ca 



The real and imaginary parts of F(p, S, g, i) remain completely unchanged and are given by Eqs. (|25| 
and (|26p . Of course, there is an important difference between the results for x and F of the last and 
the present section. In Eqs. (|52|) and (|53|) not only the coefficients but also the integration constants 
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C3, C4, C5, Cq may be arbitrary functions of q,t (similar remarks apply to F). Those conditions, from 
the set (l4"2"1) - (|51l) . which have not yet been used present constraints for the spatial and temporal variation 
of all these parameters. 

In order to find the explicit form of these constraints the results (f52"| . (fSU)) . (f2"5| . and (f2"o| have to 
be inserted in the conditions not yet used, namely (j4"4"l) , (|4"5)) (j4l))) and (|5Tj) . The calculation reported in 
detail in appendix[B]lcads to the following results for the parameters. The quantities C5, Cq are arbitrary 
functions of q, t, while C3 = C4 = 0. The complex coefficient d is an arbitrary function of q, t. The 
remaining coefficients are given by 

a=i2mud, (54) 
,5C 6 n dC b 



e=(H 1 +iH 2 )d, (56) 
where Hi is an arbitrary real function of q, t, and Hi is given by 

Bi = 2 m u(l^ + ^)-2^^-^. (57) 



<9i dt 



The quantity u, defined in Eq. (| 147|) . is an arbitrary function of t. The above relations show that all 
coefficients are proportional to d. Therefore, since C3 — C4 — 0, the coefficient d drops out of the 
differential equation. Thus, we are left with three arbitrary real functions, C5, C§, H\, of q, t and and a 
single arbitrary function u of t. All other parameters may be expressed in terms of these four functions. 

The complex quantities x an d F are given by 

x= _^ e -c V («s+c 5); (58) 
u 

F=l ^P e Ce e -<uS + C 5 )^ (59) 

md 

Inserting the results for a, b, e and dropping d the differential equation for x takes the form 

(60) 



% 2mu— (2^-12^)^ f^X 
dt \ dq dq J dq dq 2 



Hi 



, .,135, dC 6 \ n dC 5 dC 6 
Zmu — — — 



u dt dt J dq dq dq 2 



x = o 



From now on we use for simplicity the symbol for the partial derivative for all kinds of derivatives. At 
this point one may use the results (|58|) . (f59|) . to show that both sides of the fundamental condition (|4Tj) 
agree. 

In order to compare with the previous, "interaction-free" equation (I38p it is convenient to replace the 
fields u, Hi by 

p(t) = l, v- = -ct- ( 61 ) 

u Zmu 
Using these fields and rearranging terms, Eq. ([50]) takes the form 

_P_dx +2 ]^_ fdCe ^^dC 5 \ dx , p 2 d 2 x 



i dt 2m \ dq dq J dq 2m dq 2 

- v ( i dp | dc 6 \ %P * ( 2 dc 5 dc 6 { d 2 c 5 \_ Q (62) 

I \ p dt dt J 2m \ dq dq dq 2 J 

Comparing (f62l) and (|38|l one sees that now, instead of a single arbitrary function V(q,t) in (|38|) . four 
arbitrary functions p(t), V(q,t), Ce(q,t), Cs(g,i) appear in (f62|) . The field p(t) has the dimension of an 
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action and replaces H, the field V(q,t) replaces the previous mechanical potential V(q,t). Each one of 
these four functions may play, in principle, the role of an "agent" for an external interaction. However it 
is possible (and we will see in a moment that this possibility applies) that some of these functions may be 
eliminated by means of a redefinition of the variable %■ Let us try first to eliminate the functions C5, C§. 
Using the transformation 

X^X = Xe C6 " ,C5 , (63) 

Eq. (I62|) takes the form 

pdx p 2 d 2 x „_ ldp_ 

--m + 2^W~ Vx + -d-t x = > (64) 

if V is replaced by V according to the relation 

'9C fi \ 2 fdC,\ 2 d 2 C, 



v = v-l- 

2m 



dq 2 



■p-fr. (65) 



Clearly, Eq. (|d5|) is permitted because both V and V are arbitrary functions. Eq. (|64p is a generalized 
Schrodinger equation, which agrees with the standard form for time-independent p = h. Obviously, the 
functions C5, Cq are - in the framework of the present Ansatz - physically meaningless, because they may 
be eliminated from the dynamic equation by means of a simple redefinition of the state variable. 

If we use the arbitrary function V instead of V, the untransformed equation (f6"2")l takes the more 
familiar form 

p f d I dp dC 6 dC 5 , p 2 , d dC 6 dC 5 ,2 f . 

--Sm-pM + ^-V X+ ^Wq + -^- l -dp X- F * = °- ( 66 ) 

This equation is very similar to the one obtained from the standard Schrodinger equation with the help 
of a gauge transformation; the difference is (besides the time-dependent field p) that a complex function 
Cs+iCq appears instead of a real function C5. In the present context, the gauge invariance of Schrodinger's 
equation is a consequence of the inscnsitivity of the continuity equation against external disturbances. 
That part of the additional external q, t-dependence in \, that may be expressed by variable C5, Cq has 
been eliminated by the requirement that the continuity equation remains valid. As a consequence of 
this requirement new "compensating" terms appeared in the coefficients. The present derivation opens a 
slightly different view on gauge theory. This point of view seems to be new, although the central role of 
the continuity equation and its symmetries for gauge theory is of course a well-established fact. 

Can we continue this way and eliminate, as we did with C5 and Ce, the arbitrary function p(t) too from 
Eq. (|64p. i.e. replace it by a constant H ? This is apparently not the case. We may eliminate immediately 
the imaginary potential term (the one proportional to the time-derivative of p) from Eq. (|64p by means 
of the transformation 

X^Xo=XC- lnp = X -. (67) 
P 

But the result obtained this way namely 

pdxo . P 2 d 2 xo T , n , RQ , 
~^ + 2^^- VX0 = {) (68) 

does not agree with Schrodinger's equation because still each h in (f38|) is replaced by an arbitrary 
function p(t) in Eq. Of course, both equations agree if one sets approximately pit) ~ constant = h. 

Experimentally, this approximation seems to be valid; at least in the mostly investigated range of not 
too large time-inter vails. 

A time-dependence of Planck's constant H has been the subject both of theoretical speculations and 
observations over large time periods [22j . It is not clear at present wether or not such a time-dependence 
exists. If it should turn out to be real, this would rise the fascinating question of its physical origin; 
obviously it does not fit into the existing scheme to describe interactions. 
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Let us recapitulate what we have done, paying particular attention to the question of a time-dependent 
h. In the last section we chose the simplest possible Ansatz to construct Schrodinger's equation. We found 
an infinite number of equations characterized by three arbitrary constants, p, C5, Cq, and an arbitrary 
function V(q, t). Let us recall the slightly rewritten result for x, 

X = -pV/9e~ C6+iC Vf . (69) 

In the present section we tried to make room for interactions without changing the basic framework. 
Essentially the only possibility to do that (see, however, the next section) was to allow for an arbitrary 
space-time dependence of the constants p, C5 , Cg . This is analogous to the usual procedure of "gauging" , 
which means postulating a space-time dependence of the parameters of the gauge group. The standard 
gauging procedure leads then, in a next step, with the help of the concept of "compensating fields" , to 
the introduction of the gauge coupling terms. Clearly, the present approach may be considered as a 
generalized gauging procedure; the proceeding is basically the same for all three constants. One could 
argue - in favour of the reality of the time-dependent h - that, in principle, this possibility to create 
interaction may be realized by nature for all parameters if it is realized for one of the constants (we know 
that it is realized for Cq, see the next section). On the other hand, the constants C5, Cq and p in \ 
play a very different role [see Eq. (|69|) ]. Gauging of C5, Cq is compensated by additional terms in the 
coefficients, gauging of p cannot be compensated completely. Only the first p in Eq. (|69p . which plays, 
similar to Cq, the role of a scaling factor for the amplitude, may be eliminated. Combining the two above 
gauge transformations creates a state function, which is given by 

Xo = xe Ce ~ lC ^ hlp = -y/pe l i. (70) 

The variable xo still contains the field p(t). The latter presents the only real "interaction" created by the 
extended Ansatz (pf0)) . It plays, as expected, the role of an arbitrary, time-dependent scaling factor for 
the action S. 

In the rest of this work we will neglect the time-dependence of p and set p = H. Then, in the framework 
of this approximation, the extended Ansatz (|40|) of the present section did not lead to any real interactions 
and we conclude that at any rate a different, more general Ansatz is required to derive the form of the 
gauge coupling. This will be done in the next section, where all the results derived in the present section, 
will turn out to be useful. There, we will continue with the threedimensional generalization of (|66p. which 
[for p{£) — h] is given by 



h,d_ dCe_ 0Cfcv _^!_V(— d<?6 9G. 



*+^£(£ + S-«7Bx-^ = 0. (71) 



The state function x depends on g^, t and is given by 

X = -h^e- CB e<% +Cs ). (72) 

We have (for p = H) two possibilities to derive these equations. First, we expect that (f7Tj) and (f72")) may 
be derived in a straightforward manner from a three-dimensional generalization of the fundamental con- 
dition (|4ip . As a second way to derive Eq. (|71[) one may of course perform simply a gauge transformation 
of the free Schrodinger equation. We will in the rest of this work always refer to (|7T|) and (|72p as derived 
according to the first possibility, i.e. from outside quantum mechanics. The difference is important. 
Choosing the first possibility, the role of the constants C5 , C§ as "agents" of an external influence on the 
considered system is clear. Choosing the second possibility, these constants are primarily abstract group 
parameters without any direct physical significance (their role may, however, be suspected; e.g. one may 
show that C5 does not depend on the wave function [10] ) . 

5 Second attempt to derive interaction 

The plausible looking Ansatz of the last section was not successful. Trying to find a better one, it may be 
helpful to review first briefly the standard methods of introducing the (abelian) gauge field in quantum 
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mechanics. Formally, the gauge held is introduced with the help of the minimal coupling rule 

d d e r ._ . 

atirr 1 ^ (73) 

6 

where A denotes the vector potential and <f> the scalar potential. Applying ([73]) and ([74]) to Eq. ([3|) one 
obtains Schrodinger's equation for a particle in an electromagnetic held, 

h d , , , h 2 f d e 



% a^-^-^\M~ x ^ A] Vj + ^ (?5) 



If ip is again written in the form (HJ) we obtain from (|75p the two equations 



dp, d P (9S_e_ A \ = ^ (7g) 



dt dq m \ dq c 
d S 1 fdS e -A 2 K 2 A^p 

m + ^ + v =-^{W-c A ) + wT (77) 



Using the variables p and S, there is also a minimal coupling rule to create the interaction terms in (|76|) 
and (|77|) from the interaction- free equations ((5|) and (J6j) . It is obviously given by 

as dS_ _ e j. 
dq dq c 
dS dS 

How to justify the minimal coupling rule ? The standard explanation is the "principle of local gauge 
invariance" . It requires that the parameters of the global symmetry group (the constants C5 , Cq in our 
notation) be arbitrary functions of qk, t without destroying the form invariance of Schrodinger's equation. 
This requires, thus the line of argument, introduction of a compensating held, the gauge held. However, 
as Eq. (|66p shows, Schrodinger's equation is already invariant under the local group; no matter wether or 
not the new fields have any physical effect. This has been pointed out |16j a few years after publication 
of the fundamental papers [26], |24j on gauge theory. The principle of gauge invariance alone, without 
additional assumptions, cannot explain the minimal coupling rule. 

A second method to introduce a gauge field exists. There one postulates that the phase of the wave 
function (or rather a part of it) is non-integrable. This method goes back to Weyl [25] and London [12] 
and is discussed in many works. We mention only two, a compact presentation to be found in Dirac's 
monopol paper [B] and a more detailed and clearly written discussion in a book by Kaempfer [10] . The 
basic idea of the non-integrable phase may be explained very quickly by starting from Eq. (|T5|) and 
eliminating the potentials by means of a singular gauge transformation [10 . The result is an equation 

^- = -0 ^r+%WC ( 80 ) 

% ot 2m oq z 

which looks like the free Schrodinger equation but cannot be used as a differential equation because 
the wave function ipc is a multi-valued mathematical object with non-integrable phase; the notation 
ipc indicates an unspecified dependence on a path C. The gauge field may be introduced exactly the 
other way round [B]. One starts from the free Schrodinger equation, postulates the existence of a non- 
integrable phase and eliminates this singular part by means of a gauge transformation, thereby creating 
the potentials in Eq. ([75]) . How to justify the introduction of the non-integrable phase ? Apparently, this 
question has not been discussed in the literature. This is surprising since the fact that Eq. ([80]) takes the 
form of the free Schrodinger equation is very remarkable and calls for an explanation. 
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This second method to introduce a gauge held, by means of a non-integrable phase, may be adapted 
for the present problem. A possible explanation - which is however still speculative in character - for 
the non-uniqueness of the phase will also be given (in the next section). Adapting this idea means that 
we allow in the Ansatz (|40|) for a non-unique explicit space-time dependence of the state function x- 
Introducing a non-unique mathematical object in a physical theory implies more or less automatically 
a further requirement, namely that all mathematical quantities of direct physical significance are well- 
behaved (unique) functions. Such a requirement (which is for well-behaved functions obvious and not 
worth mentioning), should be formulated explicitly in the present situation. We require, in particular, 
that the state function x itself is a unique (up to regular gauge transformations) function of q, t and 
that all derivatives are unique functions of q, t; these requirements imply that we obtain a well-dehned 
differential equation with unique coefficients. Thus, our once more extended Ansatz takes the form 



X= X(P> S;q,t), 

X unique, all derivatives of x unique, 



(81) 
(82) 



where the semicolon in (|8T|) has been introduced to indicate the allowed non-uniqueness of the explicit 
(?, t-dependence. It should be pointed out that the non-unqueness is allowed but not required. It will be 
realized only if the consistency requirements (|82[) can be fulffilled (how this actually works will become 
clear in a moment). 

We have to construct, using (f8Tj) . (|82|). a state equation for x- The basic construction scheme is the 
same as in the last section; in particular we start, from a fundamental condition which agrees formally 
with fill). Thus we have to solve 



'it 



F a 



dt 



, d ,d 2 



dp dp 1 dS 
dt dq m dq 



P- 



1 d 2 S 



(83) 



together with the conditions (fSTj) and ((52")) . The problem defined by (|5T)) - ((55)) is in principle completely 
independent from the problems dealt with in the last two sections. But all the formal mathematical 
results, obtained in the last two sections, may be taken over to the present problem; the new aspect 
- the non-uniqueness - does not change these results. The non-unique q, t-dependence of x may be 
expressed by three non-unique functions, namely p{t), C$(q,t), Ce(q,t), and the derivatives of these 
functions have to be, according to condition (|82p unique functions of q, t. Thus we may immediately 
conclude that Eq. ((55)) implies relation (JHSJ) • The latter will become a well-defined differential equation, 
if we are able to transform the non-unique quantity appearing in this equation [also called x m (EH)] 
by proper manipulations into a unique function of q, t. In other words, the only remaining problem 
is to fulfill condition (f82| . Before we proceed we neglect, as discussed in the last section, in (f66| the 
time-dependence of p(t) and set p{t) = %. Further, we shall use in the following the threedimensional 
Schrodinger equation (fTT)) instead of (|66p. This is because the structures to be studied can only be 
represented adequately in a fourdimensional space-time continuum. 

Our starting point is Eq. ([72]) ; obviously, condition (|82)) is not yet fulfilled. It is also clear that the 
functions C5 and C§ may be treated separately; given that we are unable to implement uniqueness, we 
may as well eliminate one or both of these functions by means of a singular gauge transformation (or 
simply set it equal to zero). Now all relations between our functions are fixed already; there is only a 
single possibility left to implement an additional condition: The functions p and S themselfes may be 
defined to be non-unique functions of qi, t. 

First we assume that C§ cannot be implemented as a non-unique function and set Cq = in (J7TJ) and 
(|72p . The form of x shows, that the remaining non-uniqueness of C5 may be compensated by postulating 
a non-unique variable S according to 

S = -h,C 5 + hip, (84) 

where ip is a unique (up to regular gauge transformations) function of qt, t. This simple linear combination 
eliminates the non-uniqueness of C5 by means of a non-unique S in favor of a unique ip and produces a 
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unique state function x- If, instead of ip, a quantity with the dimension of an action is required, we write 
tp — S/h. Eq. (O shows that the replacement of S by S according to 

S — ► S = S + hC 5 (85) 

is compatible with the structure of the continuity equation and leads only to a redefinition of the proba- 
bility current. Otherwise, this compensation procedure would not make sense. 

Let us first follow the consequences of a non- unique C5 before we come back to Cq . In order to obtain 
an explicit representation for C5 we introduce for its derivatives, which are unique functions according 
to (|82[) . the following symbols, 

_ = ^ fe , _ = 0. (86) 
Then, C5 may be written as an integral over a path C, 

C 5 (q,t;C) = [dq' k A k (<r,t')+dt'4>(<r,t')]. (87) 

J qoM'fi 

The last formula (summation over double indices) represents one of the fundamental postulates of gauge 
theory [25], [6], [10]. The multi-valuedness of C5 is expressed by the fact that C5 does not only depend 
on the considered space-time point q, t but also on the path C which leads from a (fixed) reference 
point qo,to to q, t. All conceivable values of C5 may be obtained by specifying four real fields 4>(q,t), 
A k (q,t), k= 1,2,3. 

These four fields are not independent from each other. In the following we reproduce some well- 
known results [5], [TO] which are essential for the physical interpretation of these fields. We define a 
four-vector x^, fi — 0,1,2,3 by xq — vot, (where vq is an unknown constant with dimension of a 
velocity) and x k = ?fe, k = 1,2,3. The four-vector field is defined by Aq{x^) — (j>(qk,t)/vo and 
Ak(Xfj,) = Ak(qk, t), k = 1, 2, 3. Then, C5 may be written in the form 



C 5 (q,t;C) = C 5 (x;C) = dx'^x'). 
Integrating along an arbitrary closed path Co from x^ to x^ Stokes integral theorem implies 

where the antisymmetric tensor d/ Ml/ characterizes the infinitesimal surface element of a surface A(Cq), 
which is bounded by our curve Co and otherwise arbitrary. Given that C5 is multi- valued, the path integral 
on the left side of Eq. ([89]) must not vanish for arbitrary paths Co- This implies that the integrand of 
the surface integral, which is denoted by F^ Vl must not vanish for arbitrary space-time points. In other 
words, the set of points defined by 

must not be empty. At the points where Eq. (|90[) holds, the non-uniqueness of C5 implies the non- 
commutativity of the derivatives with regard to x^. If = for all x^, then the influence of A^ (and 
C5) may be eliminated by means of a regular gauge transformation. If, on the other hand, points exist 
where Eq. (|90p holds true, then the values F^ v takes at these points are invariant under regular gauge 
transformations. This means that F^ v has a gauge invariant meaning and implies the possibility that 
Ffj,,, plays a role as (is proportional to) a classical field or force. 

As a consequence of the definition ([90]) the field F^ v obeys the differential equation 



(91) 
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The left hand side Ta m „ of Eq. (|5Tj) is antisymmetric in all three indices. Only four relations in Eq. (f9"Tj) 
are independent. These may be written in the form 

BF 

e,x^ = 0- (92) 
dx x 

The definition of e K \^ may be found e.g. in the book by Landau [11]. If the essential components of 
F^y are renamed according to 

Fqi = —Ei, F02 — —E2, F03 = —E3, 

, (93) 

F12 = +B3, F13 = —B2, F23 = +B\ 

then the four relations l|92|) may be written as differential equations for the vektorfields E and B (with 
components Ei and Bi), 

f = „, * ; xl + ±§ = 0. (94) 

Thus, the fields E iy Bi must be solutions of the homogeneous Maxwell equations, if the constant v a is 
identified with the velocity of light c. Of course, Eqs. ([M]) are not sufficient to determine Ei, Bi; some 
more equations are required. An obvious possibility is to identify v$ with c and to postulate that the 

rest of the equations is given by the two inhomogeneous Maxwell equations, i.e. to assume that E and B 
are proportional to the electric and magnetic field vectors E und B of Maxwell's theory. This postulate, 
which is an essential part of Hermann Weyls first gauge theory, is closely related to the structure of 
Eq. ((HZD. 

We still have to find constants of proportionality with suitable dimensions. The "fields" Ei and Bi 
(both with dimension c»r 2 ) have to be connected, by means of proper constants of proportionality, to 
the two basic terms of the classical particle-field concept, force and field. If a field Ei, which has the 
dimension of a force, is defined by means of aEi = Ei, then the constant a has the dimension g 1 cm 3 see" 2 . 
Using the available constants m, h, c only a single combination with suitable dimension may be formed, 
namely he. So we set a — he. In order to have an "objective" field Ei, whose existence does not depend 
on the presence or absence of a test particle, we have to introduce one more constant, the charge e. It 
is, like the mass m, part of the description of the individual particle. Thus we write Ei = eEi, where e 
has the dimension g^ cmA sec~ x . Combining both constants, and making similar considerations for the 
magnetic field, we obtain 

~ e ~ e ~ e 

Ei = j^,Ei, Bi = ^Ei, F M „ = —F^ u . (95) 

Similar constants of proportionality occur if the potentials are replaced by the standard potentials of 
electrodynamics A^, cf>, 

~ e ~ c 

= -—^)(qk,t), A k (xp) = —A k (q k ,t). (96) 

Replacing in an analogous way C5 by C with the help of the relation 

C 5 (q,t;C) = ^C(q,t;C), (97) 

the multi-valued function C may be written as a path integral with the usual potentials appearing in the 
integrand 

C(q,t;C) = f' [dq' k A k (q',t') - ctfftfj)] . (98) 

Using these properly scaled variables, we obtain the following well-known relation between the components 
of Fpu and the potentials A k , 4>, 
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Finally, with the help of the relations 



dqk he dt h 

(and with Cq = 0) Schrodinger's equation (|102[) for the variable 



= TZ A ^ -5T = -^. ( 10 °) 



takes the form 



^-xA-^e'i S = S+-C, (101) 



"(I + + & t (l - 440V - ^ = 0. (102) 



Thus, our second attempt to derive Schrodinger's equation with gauge coupling turned out to be success- 
ful. 

Now we come back to the question if a meaningful theory may be constructed with a nonzero multi- 
valued Cg. There is a standard argument on this point in the text-book literature. It says (using the 
present notation) that a nonzero Cq is forbidden, because it leads (in contrast to a phase change) to a 
modification of the amplitude of the wave function and to a corresponding change in probability density. 
If this argument were true, a modification of the phase would be forbidden as well. The latter leads 
also to a modification of the amplitude - as a consequence of the coupling between phase and amplitude. 
However, while this argument is strictly speaking wrong it leads to the correct conclusion. We may 
proceed as we did with C5, allowing for a multi- valued p and defining a new probability density p by 
means of the relation 

p = pe- 2Ca . (103) 

The functions p and Cq are multi-valued while p is single- valued. We may represent Cq again (as we did 
with C5) as a path integral and the associated potentials would again fulfill the homogeneous Maxwell 
equations. But this is not sufficient. In order for the substitution p — ► p (together with 5* — > S) to 
make sense, it is necessary that the structure of the continuity equation remains intact. This means that 
performing the substitution p — > p in 

dp ■ 9 P f BS e -A)=0 (104) 



dt dq m \ dq c 

produces a mathematically well-defined differential equation with unique coefficients, which still has the 
structure of a continuity equation (with possibly redefined density and current). This is not the case as 
insertion of (|103|) in (|104p shows. Only the interaction mediated by a non-unique C5 can be realized in 
nature. 

6 Discussion 

In this section we will give a summary of the paper and discuss its most important results. The calculations 
performed in section[3]) led us to conjecture that Schrodinger's equation may be derived from the following 
three assumptions 

(1) The continuity equations holds for a probability density p and a probability current j, which depends 
linearly on the gradient of a function S. 

(2) The system may be described by a complex state function X- 

(3) The state function \ obeys a linear differential equation. 



We have not verified in detail this conjecture but given arguments supporting its validity in section [3]); 
in the following discussion it will be considered as true. 
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The above three assumptions define a quantization procedure. The system that has been quantized 
was, however, not a particle but rather a statistical ensemble of particles. This quantization method 
does not require a Hamiltonian; instead the conservation law of probability (first assumption) represents 
the fundamental input of this theory, which may be roughly characterized as a statistical quantization 
method. The linearity (third assumption) of the state equation may also be understood in terms of the 
statistical nature of this theory. Finally, the second assumption of a complex- valued state function is the 
essential "non-classical" part of this theory. This is a formal requirement and its physical meaning and 
origin is a priori unclear. It is, however, a very simple assumption, much simpler than e.g. the canonical 
commutation relations, which are the essential nonclassical part of the conventional quantization method 
and are also of a formal nature. 

We may ask why the present quantization method is conceptually simpler than the conventional one. 
The present method starts from a statistical ensemble while the conventional method starts from a single 
particle. Let us consider the classical limit of Schrodinger's equation. This limit may be performed 
conveniently by setting H = in Eqs. §5§ and ([6|). The result is the Hamilton- Jacobi equation {H 
and the continuity equation (|5|). The latter is obviously not eliminated by performing the limit h =>• 
because it does not contain the constant h. The classical limit (JT|), ([5j) of Schrodinger's equation is a 
field theory, which describes an infinite number of particle trajectories. Each one of these trajectories 
may be calculated for given initial conditions in the framework of classical point mechanics [there is no 
coupling to p in Eq. |T|)] , but these trajectories are only realized with a certain probability which must be 
calculated with the help of ([5]). The classical limit of Schrodinger's equation is not a particle theory but 
a statistical theory; this lends support to the ensemble interpretation of quantum theory [9]. Therefore, 
for the inverse problem, the transition from classical physics to quantum theory, a statistical law may be 
a simpler (more natural) starting point than a particle law. 

Three different functional forms for the dependence of \ on P ancl S have been used in sections [3j , SJ 
and [5]). The first, simplest Ansatz \ — xi.Pi &)■> m section [3J), has been called "interaction free". The 
result was Schrodinger's equation for a particle (ensemble) in an external mechanical potential V(q,t). 

In the second Ansatz \ — x(p, S, q, t), in section 2J, the possibility of an explicit dependence of x 
on q,t was taken into account. Each constant derived in section [3]) becomes a possible "channel" for 
this additional q, i-dependence. One of these channels leads to a new and rather exotic "influence" on 
our system, which is formally described by a time-dependent Planck constant H. All other channels are 
eliminated by additional compensating terms in the coefficients. By means of a simple redefinition of 
the state variable the resulting equation may (for neglegible variation of K) be transformed to the free 
Schrodinger equation of section[3j. This tranformation agrees exactly with the usual quantum mechanical 
gauge transformation. Therefore, the Ansatz x = q, t) contains the usual gauging procedure as a 

special case. The gauge coupling terms itself cannot be derived within this Ansatz; it turns out, however, 
that it presents nevertheless an essential step towards its introduction. 

The third Ansatz in section [5]) has been written in the form x = x{pi S] q, i), which means that a 
multi-valued external space-time dependence of x 1S permitted. The idea for this Ansatz stems from the 
well-known quantum mechanical concept of a non-integrable phase. The functional form of all depen- 
dencies is already fixed by the calculations of section [4]). It turns out that the only channel leading to 
an actual physical effect is given by the parameter C5 . Its multi-valuedness must be compensated for by 
the multi-valuedness of the quantity S (the generalization of the classical action), in order for the final 
phase of x to be well-behaved. Performing a variable substitution one obtains finally the correct gauge 
coupling terms in Schrodinger's equation. In the standard formulation of gauge theory, the compensation 
effect is meant to restore the (locally destroyed) gauge symmetry, in the present formulation it restores 
the uniqueness of the state function. The last step to the introduction of the gauge field, performed in 
section[5]), is formally nearly identical to the standard theory of non-integrable phases. There is, however, 
an important conceptual difference. The field C§{q, t) of the present theory represents by definition an 
externally controlled influence on the considered system. The multi-valuedness of the action S(q, t), which 
describes the system, is a consequence of the multi-valuedness of C^q^i). The latter may be identified 
as reason while the former is the effect of the latter. This clear cause-reason relation is absent in the 
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standard theory of the non-integrable phase. 

Let us come back to the question raised in section [5]) how to justify, or "understand" , the form of the 
gauge coupling. We know that the non- uniqueness of C5 (q, t) , which describes some external influence 
on our system, is a necessary condition for the occurrence of a gauge field. So, what is the reason for the 
non- uniqueness of C§(q, t); where does it come from ? The following is not a final answer to this question 
but rather a coherent collection of remarks intended to stimulate further research. 

The non-uniqueness of a physical quantity means that it is impossible to express its effect in a "local" 
way, by stating its value on a particular space-time point. The term "local" needs further explanation. It 
applies practically to all of the conventional scheme of formulating physical laws by means of differential 
equations; in this scheme all the information required to predict (or retrodict) the behavior of a system 
is determined by stating its values (e.g. coordinates of particles or field values at all space points) at 
a particular instant of time. There is no universal law dictating that this "local" conventional scheme 
must be always true. It may well be that instead the values of the variables in a time intervall, of finite 
extent, are required to predict the behavior of a system. Mathematically such theories are described e.g. 
by delay differential equations. 

Let us consider the classical theory of charged particles and fields, having in mind a possible breakdown 
of the conventional scheme. Nearly all applications of this theory study one of two idealized situations. 
The first is calculating the fields produced by given sources, the second is calculating the trajectories of 
particles in given fields. Both types of problems are, strictly speaking, unrealistic, even if this does not 
cause any problems in macroscopic situations. Realistic problems have to take into accout the mutual 
influence of trajectories and fields. Whenever this is done seriously [T7j, [T3], [2] one encounters systems 
of delay-differential equations which are difficult to solve and whose mathematical structure is still to 
be analyzed in detail. But if one tries to simulate a realistic "nonlocal" problem in the framework of 
the conventional scheme, one encounters necessarily problems, in the form of non-unique or multi-valued 
predictions. This is what happens already for the simplest conceivable realistic problem of this type, a 
charged particle in an external electromagnetic field under the influence of its own radiation field [19] , 

Now, quantum mechanics (as well as the present theory) uses, of course, the conventional "local" 
scheme of formulating physical laws. It is possible that the above explained "nonlocality" is - in a way 
still to be clarified in detail - responsible for the non-uniqueness of C5 (q, t) . The latter plays, according 
to the derivation reported in section [5]), a central role for the form of the gauge coupling. May we extend 
this idea about the origin of the form of the gauge coupling further, to quantum mechanics itself ? The 
above calculation was in fact a simultaneous derivation of the fundamental equation of quantum theory 
and the form of the gauge coupling. Considering it that way, we made three (or two, if we neglect the 
possibility of a time-dependent Planck constant) different derivations of quantum theory, and may ask 
which one is most fundamental. The third derivation led to Schrodinger's equation with the minimal 
coupling rule, which is obeyed by all interactions found in nature (including non-abelian generalizations 
and possibly excluding gravity). It should be considered as most fundamental, if quantum mechanics is 
to be understood as a step towards a realistic description of nature. If one adopts this point of view, 
Schrodinger's equation depends also crucially on the assumption of a non-unique external influence. We 
conclude that a detailed study of the breakdown of the "local" formulation of realistic classical physics 
may be useful in order to achieve a deeper understanding of quantum mechanics. 

This attitude towards quantum theory is not new, although a minority view. We mention, in partic- 
ular, the work of Raju [17] . He revealed clearly crucial properties of classical particle-field systems which 
have been overlooked by the scientific community for nearly a century. Of course, such findings are also 
relevant for a careful characterization of quantum mechanical nonlocality. Nelson, using a quite different 
starting point, conjectured |15j that "..quantum fluctuations may be of electromagnetic origin.." This is 
to be understood in the sense that the stochastic mechanism, which underlies quantum mechanics, could 
be due to the "non-deterministic" behavior of charged particles, as a consequence of interaction with their 
own field [15j . This conjecture is - using a different language - very similar to the above. Finally we note, 
that there is a certain overlap of the present ideas with those underlying Stochastic Electrodynamics [5] . 
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7 Concluding remarks 

In this work Schrodinger's equation with gauge coupling has been derived. The construction was based 
on the validity of a continuity equation and a linear differential equation for a complex-valued state 
variable. As an heuristic recipe this scheme may probably be extended to more general situations. 
Obvious possibilities to do that include a relativistic formulation or a multi-component version, which 
should be the analogon of non-abelian gauge theories. Other types of variables (quaternions) or even 
an application beyond the realm of special relativity seem conceivable. From a physical point of view, 
the most important generalization or modification of the present theory concerns the assumption of a 
complex state variable. This is a purely mathematical assumption which - despite of the outstanding 
structural properties of the field of complex numbers - should be replaced by an different, equivalent 
assumption which can be directly interpreted in physical terms. This is the most challenging question to 
be asked in the context of the present theory. 

A Calculation of \ and F 

Combining the derivative of (fT4| with respect to p with (fTOj) . one obtains 

ddi d\i dd 2 d\2 



dp dp dp dp 



= 0. (105) 



Using (|105|) and (fTi|) one obtains an elementary differential equation, which shows that the ratio of d 1 
and di does not depend on p; it is convenient to write this in the form 

i = e^\ (106) 
d 2 

where f(S) is an unknown function. If (|16p is multiplied by 2p and afterwards combined with (| 1 T[) . one 
obtains a differential equation 

dU 

2 P ^ = U (107) 
for the p— dependence of a quantity U(p, S), defined by 



tj( q\ d ^ d Xi dx2 , ino , 

u(p ' s) = T 2 ^s~ -as- (108) 



Inserting the solution of (|107[) in (TIT)) leads to the important intermediate result 



d 2 {p,S)=^e- fl ^, (109) 



m 



J^c^v^e-^H^. (110) 

TO 

Thus, the dependence of d\ and d 2 on p is known while the two functions h and /, which both depend 
only on S, are still to be found. 

The solution of ([TJ} and (fT5"|) with regard to and is given by 

dxi = d 2 = 8x2 = di _ , ni , 

dp d 2 d\ — dia 2 dp d 2 a\ — d\d 2 

The di may be expressed, using (fT9f and pT|) . by the df, this determines the p-dependence of the right 
hand sides of the differential equations . Integrating yields 

X2(P, g)=— 11 ^(5) +e - /(s) + H(S), (113) 
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where \d\ 2 — d\+ d\ and G(S), H(S) are unknown functions of S. In (|1 12p and (| 1 1 3|) the abbreviations 

ci = a\d\ + a 2 d 2 , c 2 = a 2 di — a\d 2l (114) 

have been used. The remaining task is the calculation of the four functions /, h, G, H. 

In order to do that we insert the solution of (|107|> for U in ()108|) and combine it with l|106|) . This 
leads to the relation 

e nS)d -M- d -§=* KS) JP- 
Using (|115[) and flSJ) we obtain the following derivatives of xi, %2 with respect to S. 

Next, we calculate the derivatives of \i un d X2 [see t|l 12[) und (|113[) ] with respect to S and equate the 
results with the derivatives as given by (|116[) und (| 1 17|) . From the resulting relation we may draw two 
conclusions. The first implies that G(S) and H(S) must be constants, say 

G(S) = C 3) = d. (118) 

The second implies two differential equations for f(S) and h(S), namely 

2m\d\ 2 d 1 _c 2 cf^+ Cl 



(119) 
(120) 



c 2 dSL x {S) oiLx{S) ' 

2ro|rf| 2 9 1 _ Cl c^ s )-c 2 

~aSL 2 (S)~ cMS) ' 
Here, the abbreviations 

Ll (S) = e~ h W (i + e 2 /( s )) , L 2 (5) = e^MS) 
have been used. After a short rearrangement, Eqs. (|119[) . (|120[) take the form 

dh df 2e 2 f \ , t \ 

-ds + asTT^-^W {c2e +C1 ^' (121) 

dh dfef-e-f 1 , _ f s ,„ . 

95 - 95eTTe^7 = 2H^ ^ Cl ) ' (122) 

which shows, that the two equations may be decoupled. Addition of Eqs. (|121[) and (|122[) yields a 
differential equation for /(<!?) alone, namely 

%'^W^ ,+ '-''>■ < 123 > 

which may be solved easily. From (|123[) and (|12ip we obtain for the functions f(S) and h(S) the final 
results 

/(S)=lntan ( \r^S + C 5 ) , (124) 



2m\d\ 

k ^=- incos + Cs ) - 4^* - ^ 

containing the constants of integration C5 , Cq and the coefficients Cj , <i; . 

This determines the functional form of the variables x an d we were looking for. Using 
Eqs. (fTT2"|) . (fTT3| . (fT24"]) . and ([T25]) one obtains (23]) and ([24]) for the real- and imaginary parts of x- For 
the real- and imaginary parts of F one obtains with the help of (f2Tjl . (|109p . and (|110p the expressions 
(1251) and (El)- 
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B Calculation of parameters 

Conditions (|44p . (143)) (|4U|) and (fSTj) have not yet been used and play the role of constraints for our 
parameters. In this appendix we insert x and F [see Eqs. (|52]l . (|53]l . (|25|) . and (J26])] in ((44]), (J45J) (|46]l 
and (|5ip and evaluate the resulting parameters. 

We introduce the following abbreviations: 

U =7^^S + C 5 , (126) 



2m\d\ 

v =4iF s + ft ' (127) 

(128) 

r=-^. (129) 

Using these abbreviations, the real and imaginary parts of x an d F are written as 

Xi=-WcosUe~ v + C 3) (130) 

X2=-TUsinf7e- y + C 4 , (131) 

F 1= Te v (di sin [/ + d 2 cos J7) , (132) 

F 2 =Te v (dicosU- d 2 smU) , (133) 

and the functions a,, 6^, d,, e,, i = 1,2 [see (fll?|) - (|22|)] take the form 

ai=Te v (d sin [/ - c 2 cost/), (134) 

^Te^ (ci cos U + c 2 sin f7) , (135) 

6i =Te v (31 sin U - g 2 cost/), (136) 

& 2 =re y (5 1 cosZ7 + g 2 sinZ7), (137) 

Ji=|d| 2 Te y sinU', (138) 

J 2 =|rf| 2 Te l/ cosC7, (139) 

e!=Te v {hisinU - h 2 cosU), (140) 

e 2 =Te v (/ii cos [/ + /i 2 sin [/), (141) 



where gi, hi are defined by 



gi = hdi + b 2 d 2 , g 2 = h 2 d x - b x d 2l (142) 
hi = eidi + e 2 d 2 , h 2 = e 2 d x - e\d 2 . (143) 



Inserting Xi, Fi in (|44[) we obtain 



ci = 0. (144) 

The latter relation occurred already in section [3]); it implies that the variable S appears only in the phase 
and not in the modulus of x- From the conditions (|45|) and (|46p we obtain the relations 

ff2+ 2|cf (s^ + ^)=0, (145) 



dq dq 

9l -2\d\ 2 ^=0, (146) 

where the abbreviation 

?1 = 4r ^ 
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has been used. Eq. (|144|) . i.e. V — C@, has already been used in (I146|) . and will also be used in the 
following to simplify the formulas. Since Eq. (|145p must hold for arbitrary S, condition (|145|) implies the 
two constraints 

d 4 = 0, (148) 



dq 



(u depends only on t) and 



dC* 



g 2 + 2\d\ 2 —^ = 



(149) 



The remaining condition (|5ip leads to a lenghty expression. We rewrite condition (|51|) in the following 
form 



T e +T a + T b + T d = 0, 

T e =eiX\ - e 2 X2, 

_ d\i _ 5x2 
T a = ai — -a 2 — , 

Tb=h— o 2 ——, 

dq dq 

- d 2 X i t d 2 X2 
Td=di—-z- - d 2 - 



(150) 
(151) 

(152) 
(153) 

(154) 



dq 2 dq 2 ' 

associating a term with each one of the coefficients e, a, 6, d. Inserting xty Fi we obtain (using ci = 0) 



T e =TWh 2 + 

Te C6 [ (C 3 /ii - C 4 h 2 ) sin U - {C A h x + C 3 h 2 ) cos U] , 
ldu dC & 
u dt dt 



T a =-T^c 2 
u 



Te Ce 

u 

Te c 6 



TT dC * TT 

c 2 -— sin (7- c 2 —- cost/ 
at dt 



9i s 



,du dC$ 

> 1 - 

dq dq 



9v 



0C3 
dq 



92- 



0C4 
dq 



-92 
sin U 



ldu dC 6 
u dq dq 



dC 3 dCA 
dq dq J 



cos U 



T d =-T^\d\ 2 
u 



U dq dq 



s du dC 5 
dq dq 



,d 2 u d 2 C 5 



dq 2 dq 2 
Te c «\d\ 2 



3 • TT dC ± TT 

- sin U — — cos U 



dq 2 



(155) 
(156) 

(157) 
(158) 
(159) 
(160) 
(161) 
(162) 
(163) 



The left hand side of Eq. (|150|) is a sum of four terms, each one beeing proportional to a linear independent 
function of p, S, namely p, pS, sin S, cos S. As a consequence Eq. (|150[) leads to maximal four sub- 
conditions. The vanishing of the coefficient of p implies 



2^-2[ — 

C2 

Jd\ 2 



C-2 



1 du dC e 
I - 

u dt dt 

dC 5 ( 1 du dC 6 



dq \ u dq 



2 


dC 5 


^ — 

C2 


9l ^q~~ 


OCA 


v d 2 C 5 l 


dq j 


' dq 2 



ldu dC e 
~ u dq dq 



= 



(164) 



The condition that the coefficient of pS vanishes does not lead to a new constraint; it is automatically 
fulfilled if Eq. (|148p holds. Finally, the condition of a vanishing coefficient of smS, cosS 1 leads to the two 
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dC 4 


( dC 3 


dC 4 


dt + 






dC 3 


( 0C 3 


dC A 


dt 




+ 9l ^a~ 



=0, (165) 



relations 

W dq 2 

-rv'i -ry_> — ! — . — ; - 'di 2 ^=() : (166) 

Now all constraints for the 12 real functions dj, bi,di, e^, C3, C4, C5, Cq are known explicitely. These 
are Eqs. |[T53)1 - |[T55|> and the above relations (fT4l)) . dUSJ), ([T35]) . ([1151) . Using (fllg)) , (fT4T)j) and (fT^5|) . 
Eq. I|164p takes the form 

a *_ a (i* + £)_** +I !4!«£_0. (167) 
c 2 \u dt dt J c 2 |a| c 2 aq^ 

With the help of ([115)1 , ([115)1 one obtains 

fca = 2m|d| u ^- + — j - 2|d| — — - ^ . (168) 

The latter formula shows that /12 is proportional to \d\ 2 , where the constant of proportionality is de- 
termined by the functions u, C5 and Cq. Exactly the same holds for ci, C2, gi, gi- The definitions of 
Ci, gi, hi show that it should be possible, for given c,-, (?;, /i,, to express the a^, bi, Ci in terms of the di. 
The resulting relations (for the complex quantities a, b, e) are given by Eqs. ([S"4")) - ([5l))) . The coefficients 
on the left hand side of Eq. (|165|) and (|166|) are arbitrary functions. This implies the (expected) trivial 
solutions 

C* 3 = 0, C 4 = (169) 

for these equations. As a consequence, the constant d drops out of the equations, the remaining arbitrary 
parameters beeing C5, Ce, Hi, u. 
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